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Abstract 



We numerically study the dynamics of an imploding hollow cylinder composed of 
dust. Since there is no cylindrical black hole in 4-dimensional spacetime with physically 
reasonable energy conditions, a collapsed dust cylinder involves a naked singularity 
accompanied by its causal future, or a fatal singularity which terminates the history of 
the whole universe. In a previous paper, the present authors have shown that if the dust 
is assumed to be composed of coUisionless particles such that these particles go through 
the symmetry axis of the cylinder, then the scalar polynomial singularity formed on the 
symmetry axis is so weak that almost all of geodesies are complete, and thus effectively 
no singularity forms by the collapse of a hollow dust cylinder. By contrast, in this 
paper, we assume that whole of the collapsed dust settles down on the symmetry axis 
by changing its equation of state. Obtained solutions are the straightforward extension 
of Morgan's null dust solution, in which no gravitational radiation is emitted. However, 
in the present case with timelike dust, infinite amount of C-energy initially stored in the 
system is released through gravitational radiation. We also show that the gravitational 
waves asymptotically behave in a self-similar manner. 
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§1. Introduction 



Gravitational collapse of massive objects and formation of spacetime singularities are 
one of the most prominent phenomena predicted by general relativity. Physical quantities 
will blow up at the spacetime singularities and thus all of known theories of physics are not 
applicable to describe physical processes realized there. This means that new physics must 
exist at the spacetime singularities. An important issue related to the attributes of the 
spacetime singularities is known as the so called cosmic censorship conjecture. There are 
two versions, the weak and strong cosmic censorship conjectures. The weak version states 
that the singularities produced by gravitational collapse are generically contained in black 
holes, whereas the strong version asserts that, generically, timelike singularities do not occur. 
A few serious counterexamples for strong version have been found, ^^"^^ while the weak version 
is still anecdotal. Thus, hereafter, we mainly focus on the weak censorship. A more precise 
formulation of the weak cosmic censorship conjecture was given by Wald.^^ 

Weak cosmic censorship conjecture: Let E he a, 3-manifold which, topologically, is the 
connected sum of TZ^ and a compact manifold. Let {hab, Kab,ip) be nonsingular, asymptoti- 
cally flat initial data on E for a solution to Einstein's equation with suitable matter (where ip 
denotes the appropriate initial data for the matter). Then, generically, the maximal Cauchy 
evolution of this data is a spacetime, {M.,gab) which is asymptotically flat at future null 
infinity with complete T"*" . 

There are two types of the counterexamples for the weak cosmic censorship conjecture; 
the first type is that effects of spacetime singularities propagating to infinities is so fatal 
that the history of the universe ceases, and the second type is that gentle physical effects 
of spacetime singularities propagate to infinities. In the first type, the causal futures of 
the spacetime singularities do not exist, and thus the spacetime will be globally hyperbolic. 
The counterexamples of this type are not the counterexamples for the strong cosmic cen- 
sorship conjecture. By contrast, in the second type, the spacetime singularities must be 
accompanied by their causal futures, i.e., domains which suffer the physical effects of the 
spacetime singularities; such singularities are called the naked singularities. This type is also 
the counterexample of the strong cosmic censorship conjecture. 

Even though we have not yet experienced fatal influences of spacetime singularities, we 
cannot assert that the first type counterexample for the weak cosmic censorship does not 
exist in our universe. As mentioned, the spacetime singularity will be a signal of the violation 
of general relativity and implies occurrence of quantum gravitational phenomena. Thus even 
if fatal influences of singularities are predicted by general relativity, the real effects can be 
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mildened by quantum effects. 

One of the most impressive example was numerically presented by Shapiro and Teukol- 
sky.^) They showed that gravitational collapse of a highly elongate mass composed of col- 
lisionless particles might be a counterexample for the weak cosmic censorship conjecture 
in accordance with the hoop conjecture which states that a black hole with horizon forms 
when and only when the mass M gets compacted into a region whose circumference in ev- 
ery direction is C < AnGM / .'^^ Exactly speaking, the numerical analyses by Shapiro and 
Teukolsky did not show that the weak cosmic censorship conjecture breaks down, but showed 
only blowing-up tendencies of physical quantities without a trapped surface. Anyway, in or- 
der to investigate whether counterexamples for the weak cosmic censorship conjecture exist, 
we need try to solve Einstein equations as a Cauchy problem to follow the time evolution 
of the 3-dimensional space like as the analyses by Shapiro and Teukolsky. However, here, 
it should be noted that we can know whether the weak cosmic censorship conjecture holds, 
only when wc specify boundary conditions on the spacctimc singularities, or in other words, 
the physical nature of the spacetime singularities. The analyses by Shapiro and Teukolsky 
lack this point. 

In this paper, we study the gravitational collapse of a hollow cylinder composed of dust. 
This system is not asymptotically fiat and thus is out of the scope of the weak cosmic 
censorship. However, here we extend the weak cosmic censorship conjcctiu^c to the spacctimes 
which have a translational invariance in one direction and the asymptotically flat nature in 
its perpendicular directions. The case of our present interest is within the scope of this 
extended weak cosmic censorship conjecture. 

No black hole (black cylinder) forms by gravitational collapse in 4-dimensional spacetime, 
if physically reasonable energy conditions liold.^^'-^*^) Thus it is widely believed that the 
gravitational collapse of a dust cylinder is a counterexample for the extended weak cosmic 
censorship conjecture. '^^^"^^^ However, it is not necessarily true. The present authors studied 
the same subject in paper I,^^^ and showed that by imposing one of physically reasonable 
boundary conditions at the symmetry axis of the cylinder, the spacetime singularity becomes 
very mild, and, as a result, almost all of geodesies are complete. The gravitational collapse of 
a hollow dust cylinder may not be a counterexample for the extended weak cosmic censorship 
conjecture. This fact will have an important meaning also for the original weak cosmic 
censorship conjecture. Wald required that the matter should be suitable such that, in any 
fixed, globally hyperbolic background spacetime (such as Minkowski spacetime), one always 
obtains globally non-singular solutions of the matter field equations starting from regular 
initial data. The dust does not satisfy Wald's requirement. However, the example of the 
collapsed hollow dust cyhnder implies that Wald's requirement is too restrictive. It is well 
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known that gravitational collapse of spherical perfect fluid can form a central shell focusing 
naked singularity-'^^)"^-'^^ which can be strong in Tipler's sense. ^^^'^^^ The perfect fluid also 
does not satisfy the Wald's requirement. However this central shell focusing singularity is 
massless,^'^) and hence it is non-trivial issue whether these examples are so serious that these 
should be regarded as counterexamples for the weak cosmic censorship conjecture. 

In paper I, we required that the dust particles do not stay on the symmetry axis, but pass 
through there. This requirement is equivalent to the assumption that the dust is composed 
of collisionless particles. By contrast, in this paper, we require that after dust particles arrive 
at the symmetry axis, these stay there. This is a straightforward extension of the Morgan's 
null dust solution, and it is one of our purposes to clarify the difference of the dust solution 
from the Morgan's null dust solution. This requirement is equivalent to the assumption that 
very high density makes interactions between particles strong. By this condition, a remnant 
of the collapsed dust cylinder remains on the symmetry axis. 

This paper is organized as follows. In §2, we show the canonical coordinate system for 
the spacetime with whole cylinder symmetry and the Einstein equations in this coordinate 
system. In §3, we briefly review the Morgan's solution which describes the collapse of a 
hollow cylinder composed of null dust. In §4, we give the basic equations for the dust 
case, and discuss the boundary conditions on the metric components and matter variables 
at the spacetime singularity. Then, in §5, numerical results are shown. In §6, we give 
analytic solutions which asymptotically agree well with the numerical solutions for the metric 
variables. This analytic solutions imply that infinite energy per unit translational Killing 
length initially stored in the system is released by gravitational radiation. Finally, §7 is 
devoted to summary and discussion. 

In this paper, we adopt the unit of c = 1. We adopt the abstract index notation; the 
Latin indices a, b, c denote the type of the tensor, while the Greek indices //., i>, p, a mean 
the components with respect to the coordinate basis. 

§2. Basic Equations for Cylindrical System 

The spacetime with whole cylinder symmetry is defined by the following metric,^^^'^^^ 

where 0<r<-|-oo,0<(^<27r and — oo < 2; < +00 constitute the cylindrical coordinate 
system, and 7, ^ and R are functions of t and r. This coordinate system is called the 
canonical coordinate. In order that r = corresponds to symmetry axis, R should vanish at 
r = 0. The coordinate variables t, r and z are all normalized so as to be dimensionless. 
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The Einstein equations for the hne element (12- II) are 
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RR' ii^^ + ij'^) - 2RRijij' + R'R" - RR' 
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- SttGv^ (^R'T\ + RT' 

RR (tp^ + - 2RR!i)^l)' + RR" - R'R' 



(2-2) 



R'^ - R? 



R " ^' 
R-R" = -SttCv^ (T*t + T%) , 

AtcG 



R 



-g {T\ + T\. + T%-T\). 



(2-3) 

(2-4) 
(2-5) 

(2-6) 



where g is the determinant of the metric tensor, and a dot represents the derivative with 
respect to t, while a dash represents the derivative with respect to r. 

§3. Imploding null dust 

It is instructive to see the gravitational collapse of an imploding hollow cylinder composed 
of null dust, before studying the case of dust. An exact solution was given by Morgan, and 
was studied by several authors.^^-''^^-''^'^) The stress-energy tensor of the null dust is given by 



-tab 



pk^k 



(3-1) 



where fc" is the vector field tangent to future-directed ingoing radial null geodesies, and p is 
assumed to be non-negative function so that physically reasonable energy conditions hold. 
Non-trivial components of geodesic equations k"-'W ak^ = are 



dkt dkr 

— — = and — — = 0, 

du du 



(3-2) 



where u is the affine parameter and we have used the null condition k°'ka = 0. For imploding 
null dust, we have 



(-1,-1,0,0). 



(3-3) 



Using the geodesic equations k"-V ak'^ = 0, the equation of motion VqT"* = becomes 

{dt-dr){Rp) = 0. (34) 

The general solution of the above equation is 

D{w) 



P 



R 



(3-5) 



6 



where w = t + r is the advanced time and D is an arbitrary non-negative function of w. 

The right hand sides of Eqs. (l2-5p and fl2-6p vanish identically for null dust. Everywhere 
finite solutions for these equations are 

,(,,.)^r^(i±£t£i^... (3.6) 



jr2 _ ^2 



and 

R = g{t + r)-g{t-r), (3-7) 

where / and g are arbitrary functions. Here, in accordance with Morgan, we assume / = 
and g{y) = y/2. Then we have 

V' = and R = r. (3-8) 

Using the above solutions, Eqs. f l2-2p and f l2-3p lead 

(dt - drh = 0, (3-9) 
{dt + drh = 87vGD. (3-10) 

From the above equations, we find that 7 is also a function of the advanced time w only, 
and we have 

7 = 47rG / D{x)dx. (3-11) 



We assume that the function D{w) has a compact support (wi, Wo) so that the null dust forms 
a hollow cylinder. Due to this assumption, 7 vanishes for w < Wi, and thus the symmetry 
axis r = is regular for t < Wi, i.e., before the null dust reaches the symmetry axis. For 
Wi < w < Wo, 7 is an increasing function of w due to the non- negative nature of D, whereas 
7 is a positive constant for w > Wq- Non- vanishing 7 at r = implies the conically singular 
symmetry axis, and thus the symmetry axis is conically singular for t > W\. 

We can easily see from Eqs. fl3-5p and fl3-8p that p is infinite at r = if D does not vanish 
there. Thus for Wi < t < Wo, the energy density p diverges at r = 0. Although all of the 
scalar polynomials composed of the Riemann tensor vanish, the components of the Riemann 
tensor with respect to a frame parallelly propagated along a timelike geodesic connected to 
the symmetry axis diverge for wi < t < Wo in the limit of r — > 0.^'^^ Hence the symmetry 
axis r = ioT Wi < t < Wo is the so-called p.p. singularity.^^ This singularity is naked and 
satisfies the strong curvature condition defined by Krolak.^''^'^^^ D vanishes on the symmetry 
axis for t > Wo, but the regularity is not recovered, since, as mentioned above, the symmetry 
axis is conically singular also for t > Wo- This conical singularity implies that a singular 
line source remains there. Thus this solution describes a process that the the null dust 
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collapses into the symmetry axis and settles down there by changing its equation of state. 
It is well known that such a singularity is obtained in the thin limit of a straight cosmic 
string, but converse is not necessarily true. As Geroch and Traschen discussed, the equation 
of state for the matter condensed into this singularity cannot be uniquely specified. 

It should be noted that all of the variables, 7, ip, R and are finite even at the spacetime 
singularity in the domain < r < +00. Thus the domain for these variables can be extended 
from 0<r<ooto— oo<r<oo. We may call the domain < r < 00 the physical domain, 
while we call the domain — 00 < r < the fictitious domain. In this extended domain, 
the condensation of the null dust into the spacetime singularity is regarded as a removal of 
the null dust from the physical domain to the fictitious one. For t < Wi, whole of the null 
dust exists in the physical domain. At time t = Wi, the inner surface of the hollow cylinder 
reaches the symmetry axis r = 0. Whole of the cylinder enters into the fictitious domain 
by t = Wo (see Fig{T]). The introduction of the fictitious domain is useful for constructing 
numerical solutions which describe the collapse of an imploding hollow cylinder composed 
of dust, in the next section. 




Fig. 1. Schematic diagram of the Morgan's null dust solution. 



It is worthwhile to note again that, as mentioned in §1, we have implicitly imposed 
some boundary condition at the naked singularity through the requirement for realizing the 
chronological future of the spacetime singularity as smooth as possible. 



8 



§4. Imploding dust 



In this section, we consider the dust whose stress-energy tensor takes the same form as 
the null dust, 

T"^ = p^/V, (4-1) 

but here -u" is the unit timelike vector field whose integration curves are the world lines of 
dust particles. Here we write the components of in the form 



u'- 



(1,^,0,0), 



and we introduce a conserved density D defined by 



D 



-gpu 



(4-2) 



(4-3) 



Note that this variable D is equivalent to D introduced in the preceding section. We also 
assume that p is non-negative so that all of the reasonable energy conditions are satisfied. 
As in the case of the null dust, we assume that D has a compact support in r-domain so 
that the dust constitutes a hollow cylinder. 

4.1. The Einstein equations 

Since, as mentioned, the metric variable it! should vanish at r = before the singularity 
formation, we write it in the form, 

R = r(3. (4-4) 

Then the Einstein equations become 

+ 2pp' + r {2p'^ + pp" - p^^ + (^p'p" - pp'^ 



+ 



p + rp' + rpv 



(4-5) 



7 = - + 2tPP' + r'(/3'' - /j^) } ' r^pp (^^^ + V^'') - 2r/3 {P + rP') ^jji/j' 
-pp + r (^pp' - pp'^ + [pp" - p'p'^ 



+ 



^P+{p + rP')v^ 



P-P"- -p' = ^e^-^Dv^T^, 
r r 



(4-6) 

(4-7) 
(4-8) 
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+ -- (l + r^^ ^' = ^e^-^nVT^. (4-9) 

jj r \ p J rjj 

The equations of motion for dust VqT"^ = lead 

D + {vD)' = Q, (4-10) 

^ = i; + W = (1 - t;2) _ + ^' _ y I . (4.11) 

Equation fl4-10p represents the mass conservation, whereas Eq. fl4-lip is the geodesic equation. 

In contrast to the null dust case, it seems to be impossible to obtain solutions for the 
above equations analytically. Thus we invoke numerical simulations to study this system. 

4.2. Boundary condition 

It is the primary purpose of this paper to numerically construct solutions for the collapse 
of the dust similar to the Morgan's solution reviewed in the preceding section. When a hollow 
cylinder composed of dust collapses to the symmetry axis, D has non-vanishing values at the 
symmetry axis. It is easily seen from Eq. fl4-3I) that the rest mass density p diverges at r = 
if D does not vanish there, as long as i? = at r = 0, i.e., /3 is finite there: as will be shown 
later, this is true in the present case. Thus, when the hollow cylinder of dust collapses to 
the symmetry axis, the spacetime singularity will form there. 

In the null dust case, (3 = 1 and if) = Q are the solutions for the Einstein equations, 
whereas these are not in the case of the dust. The motion of the dust disturbs (3 and ■?/'. Thus, 
it is non-trivial whether (3 and if) are still everywhere finite after the spacetime singularity 
appears at the symmetry axis. In the paper I, the present authors have studied the same 
subject, i.e., the collapse of a hollow cylinder composed of dust. By imposing appropriate 
boundary conditions on the metric and matter variables, /3, 7, if), D and at the spacetime 
singularity, the present authors constructed numerical solutions for these variables which are 
everywhere finite and continuous even after the formation of the spacetime singularity. At 
that time, the present authors imposed the going through boundary condition on the motion 
of the dust: all of the dust particles collapsed to the spacetime singularity again expand with 
the same speed as their collapsing speed when they reach there. 

In the present case, the same boundary condition for the metric variables, (3, 7 and ip, at 
the spacetime singularity are also available, and we can construct solutions for them. 
The chronological future of the spacetime singularity is realized and thus the resultant space- 
time singularity is naked. By contrast, for the matter variables D and f , we impose different 
boundary conditions from those assumed in the paper I. We set the boundary conditions 
on D and v so that the similar situation shown in the preceding section is realized i.e, the 
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dust collapsed to the symmetry axis settles down there; we call this boundary condition the 
settling down boundary condition. 

First, we show the boundary conditions for the metric variables at the symmetry axis 
r = 0. From Eg. (14 -8 1) , we have 

fj' = -iTrGe^f-^DVl^ + r0 - (3"). (4-12) 

The above equation gives a Neumann boundary condition on /3 at r = as 



(3' = -ArrGe'^-^'DVT^ l^^^. (4-13) 



By the same procedure, we obtain the Neumann boundary conditions on and 7 at r = 
from Eqs. and (gS]): 



7 ■ 



47rG 



(4-14) 



^' = — ^e^-^Dv/r^ (4-15) 

where we have used Eq.f l4-13l) to derive Eq. fl4-14l) . These boundary conditions guarantee 
the finiteness of /?, 7, and their derivatives with respect to r as long as D is finite and 
f ^ < 1. In other words, the spacetime singularity formed at the symmetry axis r = causes 
at most regular singular points in the Einstein equations as long as D is finite and < 1. 
If 7, 7', 'i/' and 'i/'' do not diverge, the right hand side of Eq. fl4-lip is also finite. Therefore, 
if initial data is appropriate, D will remain finite and will be always less than unity even 
at the spacetime singularity. As a result, we will have finite and continuous solutions for all 
variables /3, 7, V', D and v even at the spacetime singularity. 

Here we should note that non-vanishing derivatives of the metric variables with respect 
to r at the symmetry axis imply the irregularity of the spacetime at the symmetry axis, even 
if the metric variables and thus the components of the metric tensor are functions of t 
and r. The functional regularity of the metric components is not equivalent to the spacetime 
reg ularity.^^) It is seen from Eqs.( l4-13|l - fl4-14p that if D does not vanish at the symmetry 
axis r = 0, the symmetry axis becomes singular. 

As mentioned, in order to determine the boundary conditions on the matter variables D 
and V at the spacetime singularity, we refer the Morgan's null dust solution shown in §3. 
We extend the domain < r < 00 to —00 < r < 00; as in the case of the null dust, we call 
the original domain the physical domain and the additional domain the fictitious domain. 
We assume appropriate 7 and ip in this fictitious domain r < 0; as in the paper I, we may 
assume 7(t, r) = 7(t, — r) and ip{t, r) = ip{t, — r). Then we solve Eqs.f l4-10|) and (14-lip for the 
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Fig. 2. 



Initial Hypersurface 

Schematic diagram of the numerical domain. 



fictitious domain as well as for the physical domain. We set the initial condition so that all 
of the collapsing dust will enter into the fictitious domain from the physical domain. When 
the dust passes through r = 0, the symmetry axis becomes the spacetime singularity since 
p and thus the Ricci scalar diverge there. In contrast to the null dust case, it is nontrivial 
whether there may remain the (^-angular deficit, after whole of the dust enters into the 
fictitious domain, or in other words, after whole of the dust is condensed into the spacetime 
singularity. By investigating Eg. (14-61) . we find that the conically singular symmetry axis 
remains as a final product by the collapse of the hollow cylinder composed of dust. Equation 
fl4-6l) leads an evolution equation for 7 at r = as 



(4-16) 



Note that P'^e"^ should be unity on the regular axis. However, from the above equation, 
we can see that even if (3~^e'^ is initially unity, it becomes larger than unity after the dust 
collapses to the symmetry axis r = 0, since f < for collapsing dust. (To recognize the 
situation, see FigEl) 

§5. Numerical Simulations 

5.1. Initial Data and C-energy 

We set the initial conserved density and velocity field as 

15a 



D 



327rw5/5 



r -l{l-w)f[r -l{l + w)f . 



(5-1) 
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l-exp(-^), (5-2) 

for l{l — w) < r < l{l + w) and vanishes elsewhere, where a is the rest mass per unit KiUing 
length in the direction with translational invariance, / is a positive parameter to specify the 
mean-radius of the hollow dust cylinder, w is a positive parameter smaller than unity, which 
specifies the thickness of the hollow dust cylinder, and fi is the parameter to control the 
gradient of the velocity field v. 

We set the initial data of metric variables, /3, 7 and and their time derivatives in the 
following manner. We set /5 = 1 and (3 = 0. Then the constraint equations (14- 5 p and (14- 6 p 
become 

1 =rip + — , (5-3) 

V 1 — v'^ 

7 = r. 2 • (^'^^ 



In order to obtain 7, we have got to integrate numerically Eg. (15-31) . whereas Eq.f l5-4p gives 
directly the time derivative of 7. We are interested in the initial situation similar to the static 
configuration as well as possible although the initial imploding velocity does not vanish. Thus 
we set tp = In order to determine the initial data of ■?/', we use Eq. fl4-9p with ^Z' = 0, 

^P" = -1 {^^' + AnGe''-'^ dVi^^ ■ (5-5) 

We numerically integrate Eqs. (15-30 and (15-50 simultaneously outward from r = by imposing 
the boundary conditions ip\r=o = = 4''\r=o and 7|r=o = which guarantee the regularity of 
the initial data. 

The vacuum region of the initial data obtained by the above procedure agrees with the 
Levi-Civita solution, 

ip = —nlnr, 7 = fi;^lnr + A and = 1, (5'6) 

where k, A are constant numbers which characterize this solution. Integrating Eq. (15-51) . we 
find that k vanishes for r < 1{1 — w) . Then from the regularity condition 7|r-=o = 0, we have 
A = for < r < /(I — w). For r > /(I + w), we have 

rlil+w) 

K = 47cG / dre^-'^DVl - v^. (5-7) 

Jl{l-w) 

Since D is positive for 1{1 — w) < r < Z(l + w), k is positive in the domain r > Z(l + w). 
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The C-energy Eo{t,r) proposed by Thome is the quasi-local energy per unit Killing 
length in the direction with the translational invariance for the spacetime with whole cylinder 
symmetry.^^^ Its definition is given by 



(5- 



Substituting Eq. (l5-6p into the above equation, we have 



Eo{t = 0,r) = ^{KHnr + X) (5-9) 

for the vacuum region, < r < /(I —w) or r > l{l+w), of the present initial data. Since both 
of K, and A vanish in the inside vacuum region, the C-energy vanishes for < r < Z(l — w). 
The total C-energy is obtained by taking the limit of r — > oo. We can easily see that the 
total C-energy of the present initial data is infinite irrespective of the values of A and k. 
This means that the total C-energy is infinite irrespective of D and v, but this fact is not 
so terrible. The similar situation is realized also in the Newtonian cylindrically symmetric 
system; the logarithmic divergence of the Newtonian gravitational potential at the spatial 
infinity necessarily leads infinite gravitational binding energy. 

In order to make the total energy per unit translational Killing length finite for the 
Levi-Civita spacetime, Thorne also proposed another definition of the C-energy as 

i^^N = ^ (1 - e-«^^o) = ^ [l + e-- {{r$r - {r^f}] . (5-10) 

Substituting Eq. (l5-6p into the above equation, we have 

i^N = ^ (l - e-V--^) (541) 

By taking a limit of r ^ oo in the above equation, we find that the total value of is equal 
to 1/8G irrespective of the values of A and k. As shown by Hayward, 1/8C is the upper 
bound of total value of the new C-energy if the null energy condition is satisfied and if there 
is no singularity in the initial data.^^-* 

5.2. Evolution 

In order to study the dynamical behavior of the dust and spacetime geometry, we numer- 
ically integrate Eqs. fl4-7p - fl4-lip . We adopt the finite difference method and MacCormack 
scheme to solve the equations for the metric variables Eqs. (14- 7p - (14-91) . In order to solve the 
equations of motion for the dust (14-lOp and (14- lip , we adopt the method invented by Shapiro 
and Teukolsky; we follow the motion of large numbers of cylindrical mass shells which move 
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along timelike geodesies and then construct the conserved rest mass density D and the ve- 
locity field V from their positions and velocities.*-' We show numerical solutions for P, ip, 7, 
D and Dv with the parameters I = 1, w = 0.5, a = 10~^ and fi = 10~^ in FigsJSHTl The 
numerically covered domain is < r < 15. The numbers of the spatial grid points and the 
mass shells are 3 x 10^ and 300, respectively. 

The constraint equations fl4-5l) and fl4-6l) are satisfied by the numerical solutions for 
Eqs. (14- 7p - (14-91) only if these numerical solutions are good approximation of true solutions 
for the Einstein equations. Thus Eqs. (14-51) and (14 -6p are used for monitoring the accuracy 
of numerical integrations of Eqs. (14- 7p - (14-91) . The relative errors estimated by the constraint 
equations are less than 10"'^ in the numerical data shown in these figures. 

1.06 
1.05 
1.04 
1.03 
1.02 
1.01 

1 

2 4 6 8 10 12 14 16 

r 

Fig. 3. Several snapshots of /3 in the physical domain r > 0. The horizontal axis represents the 
radial coordinate r. 

The inner surface of the hollow cylinder reaches the symmetry axis r = at t ~ 3.0. 
As mentioned, when the dust reaches the symmetry axis r = 0, the rest mass density p 
and thus the Ricci scalar blow up there, since (3 is finite there. This implies that the s.p. 
curvature singularity^^ forms there. Since, as can be seen from these figures, all of the metric 
variables /3, 7 and ip are everywhere finite and continuous, there is the chronological future 
of this singularity. Therefore a naked singularity forms in the spacetime constructed by this 
numerical simulation. 

§6. Asymptotic Behavior and Physical Implication 

The numerical simulations showed that a ripple in the metric variable (3 propagates to 
infinity and /5 asymptotically approaches to unity. By contrast, ip and 7 show different 
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Fig. 4. The same as Figj3l but for ■0- 
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Fig. 5. The same as FigjSl but for 7. 

behaviors. After whole of the dust cyhnder is condensed into the spacetime singularity, any 
characteristic scales disappear in this system. Thus the self-similar behaviors are expected 
for the metric variables ip and 7. In order to obtain asymptotic solutions for if) and 7 
analytically, we introduce a variable defined by 



t - L 



(6-1) 



where tg is a constant. Then the metric variables ip and 7 are expected to asymptotically 
depend on only ^, but, as shown below, this is not true. Assuming that P = 1 and ip depends 
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Fig. 6. Several snapshots of the conserved rest mass density D in the physical domain. The 
horizontal axis also represents the radial coordinate r. Note that the depicted range is different 
from those of FigsiSHS] 
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Fig. 7. The same as FiglH but for the momentum density Dv. 



on only ^, Eq. fl4-9p becomes 



(ftp dtp 



If we also assume that 7 depends on only ^, we have, from Eqs. ( 14-51 ) and ( 14-6p . 



(6-2) 



dip\' 



dC 



(6-3) 
(6-4) 
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The above two equations lead dip / dC, = = d'y / dC,; there are no non-trivial vacuum solutions 
for ip and 7 which depend on only ^. 

We construct asymptotic solutions in the following manner. First, we solve Eg. (16-21) . We 
have for |^| < 1 

i/j = -KsSin'^^ + ips, (6-5) 

and for |,^| > 1 



ip = —Ks In 



(6-6) 



where Kg and ips are integration constants. Since we are interested in the late time asymptotic 
behavior, the solution (16-51) is inappropriate. The solution (16-61) diverges logarithmically at 
r = and thus this solution itself also is not what we need here. However, since the evolution 
equation for ip is linear, the solution which is finite at r = is obtained by superposing the 
Levi-Civita solution on the solution (16-61) : 



-Ks In 



In r + ips 



-Kg In 



+ ^s. (6-7) 



Next we assume that 7 depends on only ^. Substituting Eg. (16 -71) into Egs. (l4-5p and (14-61) . 
we have an identical eguation 

rf7_ 2^s(v^^~0 
d^ 



e-1 



(6- 



We can easily integrate the above eguation and obtain 

7 = 2kI In 



2Ve 



1 



+ As 



(6-9) 



where As is an integration constant. Egs. (16-71) and (16-91) are the solutions that we need. The 
Kretschimann invariant of this spacetime is given by 



X 



(1 - Ks - 2kI) r + (2 + Ks + 2k2) 



(6-10) 



where t = 1 — 1^. At r = |t — ts|, the Kretschimann invariant diverges in the case of k^ < 3/8. 
All of the numerical solutions presented in this paper satisfy k^ < 3/8. Rigorous derivation 
of this solution from a view point of self-similarity is given by two of the present authors, 
TH and KN, and their collaborator Nolan. 

The solutions (16 -7p and (I6-9P are depicted in FigsJH] and [9] together with numerical solu- 
tions; the numerically covered domain is < r < 30, and the numbers of the grid points 
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Fig. 8. Several snapshots of the metric variable ip are depicted; white circles represent the numer- 
ical values, whereas the solid curves represent the analytic solution ()6-7p . Note the plotted 
range. 
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Fig. 9. The same as FiglHl but for 7. The plotted range for each snapshot is restricted to the 
neighborhood of the symmetry axis r = 0. 

and the mass shells are 6 x 10'^ and 300, respectively. We set the parameters and As to be 
equal to the numerical values k = 2.06 x 10~^ and 7|t=3o = 4.21 x 10~^, respectively. Then 
we set ts = 5.30 so that the analytic solutions agree well with numerical data. The analytic 
solution is available only for < r < \t — ts\, and hence we have plotted the data for this 
domain. It is seen from these figures that the numerical solutions asymptotically approach 
to these analytic solutions in the neighborhood of the symmetry axis r = 0. Therefore, even 
if we do not invoke long time numerical simulations, we can know the asymptotic behavior 
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through these analytic solutions; ip monotonically decreases, whereas 7 approaches to Ag. 
Since the ripples in f3 propagate outward in the manner /? ~ 1 + f{t — r)/r in late time, 
where f{x) is a function of compact support, the C-energy at any finite radial coordinate r 
has a following limit 

limEN = ^(l-e-2^=). (6-11) 

Since the final E-^ is constant, the C-energy concentrates to the symmetry axis in the final 
configuration. 

The distance from the symmetry axis to a point labeled by a non-vanishing radial coor- 
dinate r becomes larger as time goes on, because ip —00 and 7 — As for t ^ 00. Further 
the Riemann tensor R^'^pa behaves as ^-^{i+ks) asymptotically at any radial coordinate r 
(see Appendix). Thus the final spacetime is flat except at the symmetry axis r = which is 
conically singular. Due to the settling down boundary condition, the remnant of the collapse 
of an imploding hollow cylinder of dust is the same as that of the null dust. 

The total energy per unit translational Killing length decreases by the emission of gravi- 
tational radiation. The total value of the new C-energy is initially equal to 1/8G, while it 
finally becomes (1 — e~^^=)/8G. The energy e~^'*'"/8G has been released by the gravitational 
radiation. The numerical results imply that As ~ Aa/G for 10~^ ^ o-fG < 10~^ (see Fig JTUl) . 
Accordingly the ratio of the emitted energy to the initial one e = e'"^^^ depends on a in the 
manner 

e - e-^"l^. (6-12) 

The lighter the weight of the dust cylinder is, the larger the efficiency e becomes. In zero- 
mass limit, the efficiency e becomes unity. This seems to be paradoxical, but we should note 
that it is a non-trivial matter how to relate the present results with the asymptotically fiat 
cases, e.g., the situation treated by Shapiro and Teukolsky.*^) From the point of view of the 
original C-energy Eq^ the infinite amount of energy is released by the collapse of the dust 
cylinder; we can easily see that total value of the original C-energy finally becomes 

hm Eo = ^ (6-13) 

for any radial coordinate r, whereas it is initially infinite. 

§7. Summary and Discussion 

We constructed numerical solutions for the Einstein equations, which describe the collapse 
of an imploding hollow cylinder composed of dust with a requirement that the dust particles 
stay the symmetry axis after these reach there. A spacetime singularity forms at its symmetry 
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Fig. 10. The relation between Ag and a is depicted. 



axis. Although the rest mass density and the curvature polynomials blow up at the spacetime 
singularity, this spacetime singularity causes at most regular singular points in Einstein 
equations. Thus if appropriate boundary conditions are imposed, components of the metric 
tensor are everywhere finite. Then the causal future of the spacetime singularity exists and 
the resultant spacetime singularity is naked. 

We also obtained an analytic solution which asymptotically well agrees with our nu- 
merical solutions. This asymptotic solution reveals that infinite amount of energy per unit 
translational Killing length is released to infinity by the gravitational radiation, and a con- 
ical singularity remains at the symmetry final product. Strictly speaking, this is 
a counterexample of the second type for the weak cosmic censorship conjecture. However, 
since this naked singularity is merely conical, it is not so serious. 

It might be a surprising fact for some readers that the remnant naked singularity formed 
by gravitational collapse of a dust cylinder is weak. In Newtonian theory of gravity, collapse 
of a dust cylinder finally produces a gravitational potential which logarithmically diverges at 
its symmetry axis, and thus the singularity is strong in Tipler's sense. The reason of this dif- 
ference between Newtonian gravity and relativity is as follows. The Newtonian gravitational 
potential ^ produced by a dust cylinder of radius r = £ is 




for r > £, 



(7-1) 



where Tc is an integration constant, and a is the mass per unit length, 




(7-2) 
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Usually, we assume that a is conserved in the framework of Newtonian gravity. Then ^ 
logarithmically diverges at r = when the dust cylinder becomes infinitesimally thin, i.e., I = 

0. By contrast, the present prescription does not guarantee the constancy of o", and require 
that a vanishes finally. Thus, if the present prescription is adopted in the framework of 
Newtonian gravity, the Newtonian gravitational potential finally vanishes. This prescription 
seems to be unphysical from a point of view of the mass conservation in Newtonian theory. 
However, it is not so from a point of view of the C-energy conservation in general relativity. 
Although the tidal force finally vanishes except on the symmetry axis, the C-energy is 
condensed on the symmetry axis, and produces a conical singularity there, in the framework 
of general relativity. The present prescription requires that the equation of state changes 
from dust to something other than dust, and the Newtonian approximation is not applicable 
to this changed equation of state. 

Here it might be useful to compare the present analysis with that of paper I. In paper 

1, we required that the dust particles pass through the symmetry axis after these reach 
there. In the resultant spacetime, almost all of geodesies are complete and thus this may not 
be a counterexample for the extended weak cosmic censorship conjecture, as mentioned in 
Sec.l. This example revealed that the gravity produced by the collapsed dust cylinder is too 
weak to confine the collisionless particles to the symmetry axis. In contrast, in the present 
prescription, matter is confined to the symmetry axis, but it is due to not the gravity, but the 
change of interactions between particles, or the change of the equation of state. It should 
be noted that this change in the equation of state is a result by requiring the smoothest 
behaviors of the metric components like as the Morgan's null dust solution. 

Finally, we note that, in order to complete an analysis to see whether a numerical model 
like as that of Shapiro and Teukolsky is a counterexample for the cosmic censorship con- 
jecture, we need to specify the boundary conditions at the "spacetime singularity". If we 
do not so, we can not know whether the "singularity" is really singular, or whether the 
cosmic censorship conjecture holds in a real singularity case. For this purpose, we need the 
knowledge of the global structure of the spacetime. In this sense, the analyses of paper I 
and the present paper are the first step of the attempt toward the numerical study of the 
weak cosmic censorship conjecture. 
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Appendix A 

Riemann tensor of self-similar gravitational waves 



The components of the Riemann tensor of the self-similar gravitational waves (16-71) and 
fl6-9p are given as follows: 



-(4k2+2ks+1) 



-(4k2+2k3+1) 



(A-1) 



X 



:i - K,y + (1 + k,)tVt^ - r2 - 



tzi 



(A-2) 



-(4k2+2ks+1) 



X 



R 



tz 



X 



Kst'^ - (1 + Ks)TVr2 - - (1 + /ts)?' 

-2^"'/ts(l + 2K,)e2(^=-^=)r(r2 - r2)2-=-3/2 (^^ ^ ^^^^72 
(l-2/€3)r + (l + /ts)Vr2-r2 



(A-3) 



-(4k2+2ks+1) 
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